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Abstract 

We study £)+ decays to final states involving the 77 ' with a 482 pb“^ data sample collected at ^/s = 4.009 GeV 
with the BESIII detector at the BEPCII collider. We measure the branching fractions B{Df —>■ t]'X) = 
(8.8±1.8±0.5)% and B{Df —>• = (5.8 ± 1.4 ± 0.4)% where the first uncertainty is statistical and 

the second is systematic. In addition, we estimate an upper limit on the non-resonant branching ratio 
BiDt ^ ? 7 ' 7 r+ 7 r°) < 5.1% at the 90% confidence level. Our results are consistent with CLEO’s recent 
measurements and help to resolve the disagreement between the theoretical prediction and CLEO’s previous 
measurement of B{Df —>■ 

Keywords: BESIII, Ds, Branching Fractions 


1. Introduction 

Hadronic weak decays of charmed mesons provide important information on flavor mixing, CP violation, 
and strong-interaction effects [1]. There are several proposed QCD-derived theoretical approaches to handle 
heavy meson decays [2, 3, 4, 5, 6]. However, in contrast to B mesons, theoretical treatment of charmed 
mesons suffers from large uncertainties since the c quark mass is too light for good convergence of the 
heavy quark expansion but still much too massive for chiral perturbative theory to be applicable. Currently, 
theoretical results for the partial decay widths of ground-state charmed mesons agree fairly well with ex¬ 
perimental results. However, there exists a contradiction concerning the branching fraction B{D'^ —>■ p'p'^). 
CLEO reported (12.5 ± 2.2)% [7], while a generalized factorization method [8] predicts a factor of four less, 
(3.0 ± 0.5)%. Summing the large experimental value of B{Df —>■ rj'p^) with other exclusive rates involving 
r]' gives B{Df —>• rj'X) = (18.6 ± 2.3)% [9], while the measured inclusive decay rate B{Df —>• p'X) is 
much lower, (11.7 ± 1.8)% [10], where X denotes all possible combinations of states. Therefore, further 
experimental study of the p' decay modes is of great importance for resolving this conflict. 

Recently, CLEO reported an updated measurement of B{Df —>• 77'7r+7r°) = (5.6 ± 0.5 ± 0.6)% [11]; this 
includes the resonant process p'p'^. This is much smaller than the previous result [7]. In this paper, we 
report the measurements of the inclusive rate B{Df p'X) and the exclusive rate B{Df —>• p'p'^) at the 
BESIII experiment. 

2. Data Sample And Detector 

The analysis is carried out using a sample of 482 pb“^ [12] e+e” collision data collected with the BESIII 
detector at the center of mass energy i/s = 4.009 GeV. 


Preprint submitted to Physics Letters B 


October 13, 2015 




The BESIII detector, as described in detail in Ref. [13], has a geometrical acceptance of 93% of the solid 
angle. A small-cell helium-based main drift chamber (MDC) immersed in a 1T magnetic field measures 
the momentum of charged particles with a resolution of 0.5% at 1 GeV/c. The electromagnetic calorimeter 
(EMC) detects photons with a resolution of 2.5% (5%) at an energy of 1 GeV in the barrel (end cap) region. 
A time-of-flight system (TOE) assists in particle identification (PID) with a time resolution of 80 ps (110 ps) 
in the barrel (end cap) region. Our PID methods combine the TOE information with the specific energy 
loss (dE/dx) measurements of charged particles in the MDC to form a likelihood £{h){h = tt,K) for each 
hadron (h) hypothesis. 

A GEANTd-based [14] Monte Carlo (MC) simulation software, which includes the geometric description of 
the BESIII detector and the detector response, is used to optimize the event selection criteria, determine the 
detection efficiency and estimate background contributions. The simulation includes the beam energy spread 
and initial-state radiation (ISR), implemented with kkmc [15]. Allowing for a maximum ISR photon energy 
of 72MeV, open charm processes are simulated from D'^Dj threshold at 3.937 GeV to the center-of-mass 
energy 4.009 GeV. Cross sections have been taken from Ref. [16]. For background contribution studies and 
the validation of the analysis procedure, an inclusive MC sample corresponding to an integrated luminosity 
of 10 fb“^ is analyzed. In addition to the open charm modes, this sample includes ISR production, continuum 
light quark production and QED events. The known decay modes are generated with evtgen [17] with 
branching fractions set to the world average values [9], and the remaining unknown events are generated 
with LUNDCHARM [18]. 


3. Data Analysis 

3.1. Measurement of B{D^ —>• t]'X) 

For data taken at 4.009 GeV, energy conservation prohibits any additional hadrons accompanying the 
production of a Df Dj pair. Following a technique first introduced by the MARK III Collaboration [19], the 
inclusive decay rate of Df —^ r]'X is measured. We select single tag (ST) events in which at least one 11+ or 
Dj candidate is reconstructed, and double tag (DT) events in which both Df and Dj are reconstructed. To 
illustrate the method, we take the ST mode D~ —>■ a and the signal mode Df —>■ rj'X for example. The rj' 
candidates in the signal mode are reconstructed from the decay mode r]' —> 7r+7r“?7 with the ry subsequently 
decaying into 77 . The ST yields are given as 


ysT = ^D+DT^iDs ^a)egT, 


( 1 ) 


where is the number of produced pairs and eg^p is the detection efficiency of reconstructing 

D~ —>• a. Similarly, the DT yields are given as 


2 /St = ^ a)B{Dt ^ 11 'X)B^^P^ef,^, ( 2 ) 

where B^P^ is the product branching fractions B{r]' —>■ Tr+Tr^ry) • B{ri —>■ 77 ), egrp is the detection efficiency 
of reconstructing Dj —>■ a and Df —>■ rj'X at the same time. With egrp and egrp estimated from MC 
simulations, the ratio of y/grp to t/grp provides a measurement of B{Df —>• jy'A), 


B{Dt ^ PX)B 


PDG 

7 )' 


ysT 


^ST 

^DT 


( 3 ) 


When multiple ST modes are used, the branching fraction is determined as 


B{Dt ^ PX)BlP^ = 


2/dT 
Sa 2 /sT 


2/dt 


=-nT 

^ST 


2/ST 


=-DT 

£gT 


( 4 ) 


where j/dt = Ya 2/dt i® total number of DT events. 

In this analysis, the ST events are selected by reconstructing a Dj in nine different decay modes: KgK~, 

K~Tr~, Ar+iF“7r“7r®, KgK'^Tr~Tr~, 7r+7r“7r“, 7r“ry, 'K~rj'{rj' —^ 7r+7r“ry), 7r“ry'(7y' —^ p®7,p® —^ 7r+7r“), 
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Table 1: Requirements on AE for ST Ds candidates. 


ST mode a 

data (GeV) 

MC (GeV) 

K^sK- 

(-0.027,0.021) 

(-0.025,0.021) 

K+K-n- 

(-0.032,0.023) 

(-0.031,0.024) 

K~n~TT^ 

(-0.041,0.022) 

(-0.041,0.022) 

KgK~^TT~TT~ 

(-0.035,0.024) 

(-0.032,0.026) 

TT'^'K~TT~ 

(-0.036,0.023) 

(-0.033,0.025) 

7T~r] 

(-0.038,0.037) 

(-0.041,0.032) 

VlT-TTri 

(-0.035,0.027) 

(-0.034,0.028) 

■^~Vp-y 

(-0.035,0.022) 

(-0.035,0.021) 

'!T~Tr^r] 

(-0.053,0.030) 

(-0.053,0.028) 


and T:~'K^rj. The DT events are selected by further reconstructing an ry' among the remaining particles not 
used in the ST reconstruction. Throughout the paper, charged-conjugate modes are always implied. 

For each charged track (except for those used for reconstructing K'g decays), the polar angle in the MDC 
must satisfy |cos0| < 0.93, and the point of closest approach to the e+e“ interaction point (IP) must be 
within ±10 cm along the beam direction and within 1 cm in the plane perpendicular to the beam direction. 
A charged K{tt) meson is identified by requiring the PID likelihood to satisfy C{K) > C{tt) {C{'k) > C{K)). 

Showers identified as photon candidates must satisfy the following requirements. The deposited energy 
in the EMC is required to be larger than 25MeV in the barrel region (| cos0| < 0.8) or larger than 50MeV 
in the end cap region (0.86 < | cos0| < 0.92). To suppress electronic noise and energy deposits unrelated 
to the event, the EMC time deviation from the event start time is required to be 0 < T < 700 ns. Photon 
candidates must be separated by at least 10 degrees from the extrapolated positions of any charged tracks 
in the EMC. 

The ATg candidates are formed from pairs of oppositely charged tracks. For these two tracks, the polar 
angles in the MDC must satisfy | cos0| < 0.93, and the point of closest approach to the IP must be within 
±20 cm along the beam direction. No requirements on the distance of closest approach in the transverse 
plane or on particle identification criteria are applied to the tracks. Their invariant mass is required to 
satisfy 0.487 < < 0.511 GeV/c^. The two tracks are constrained to originate from a common 

decay vertex, which is required to be separated from the IP by a decay length of at least twice the vertex 
resolution. 

The TT*^ and ?7 candidates are reconstructed from photon pairs. The invariant mass is required to satisfy 
0.115 < M( 77 ) < 0.150 GeV/c^ for 7 r°, and 0.510 < < 0.570 GeV/c^ for 77. To improve the mass 

resolution, a mass-constrained fit to the nominal mass of tt*^ or ij [9] is applied to the photon pairs. For 
77' candidates, the invariant mass must satisfy 0.943 < < 0.973 GeV/c^ and 0.932 < < 

0.980 GeV/c^. For the 77 ),.^ candidates, we additionally require 0.570 < M( 7 r+ 7 r“) < 0.970 GeV/c^ to reduce 
contributions from combinatorial background. 

We define the energy difference, AE = E — Eq, where E is the total measured energy of the particles 
in the Dj candidate and Eq is the beam energy. The D~ candidates are rejected if they fail to pass AE 
requirements corresponding to 3 times the resolution, as given in Table 1. To reduce systematic uncertainty, 
we apply different requirements on AE for data and MC samples. If there is more than one Dj candidate 
in a specific ST mode, the candidate with the smallest |Ai5| is kept for further analysis. 

To identify ST signals, the beam-constrained mass Mbc is used. This is the mass of the Dj candidate 
calculated by substituting the beam energy Eq for the measured energy of the Dj candidate: = 

Eq — where p is the measured momentum of the Dj candidate. True Dj —>■ a single-tags peak at the 
nominal Dj mass in Mbc- 

We fit the Mbc distribution of each mode a to obtain j/grp. Background contributions for each mode are 
well described by the ARGUS function [20], as verified with MC simulations. The signal distributions are 
modeled by a MC-derived signal shape convoluted with a Gaussian function whose parameters are left free 
in the fit. The Gaussian function compensates the resolution difference between data and MC simulation. 
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Figure 1: Fits to the Mbc distributions for the ST Ds candidates. In each plot, the points with error bars are data, the dashed 
curve is the background contribution and the solid line shows the total fit. 


Figure 1 shows the fits to the Mbc distributions in data; the fitted ST yields are presented in Table 2 along 
with the detection efficiencies estimated based on MC simulations. 

To select events where the decays to rj'X, we require that the DT events contain an rj' candidate 
among the particles recoiling against the ST candidate. As mentioned above, the r]' candidates are recon¬ 
structed in the decay rj' —>■ with the rj subsequently decaying into 77 . All particles used in the ij' 

reconstruction must satisfy the requirements detailed above. If there is more than one rj' candidate, the 
one with the smallest AM = \M{rj'^^^) — rn{rj')\ is kept, where rn{rj') is the nominal rj' mass [9]. The decay 
mode rj' —>■ p °7 is not used due to large contributions from combinatorial background. 

There are peaking background contributions in Miij'^^^) produced by events in which there is a wrongly- 
reconstructed D~ tag accompanied by a real rj' in the rest of the event. To obtain the DT yields, we 
therefore perform a two-dimensional unbinned fit to the variables Mbc(q^) and M{rj'^^^). For MBc(<a), the 
fit functions are the same as those used in the extraction of j/grp. For M{rj'^^^)^ the signal is described 
by the convolution of a MC-derived signal shape and a Gaussian function with parameters left free in the 
fit. Background contributions in Mijj'^^^) consist of (a) D+Dj events in which Dj decays to the desired 
ST modes, but the D+ decay does not involve an 7 '; {h) other (non-ST signal) decays of D~ and also 
non-D^ D~ processes. Component (a) is described with a first-order polynomial function. Component 
(&) is modeled with the sum of two Gaussian functions plus a quadratic polynomial function. The means 
of the two Gaussians are fixed to the rj' nominal mass [9]. Other parameters and all the amplitudes are 
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Figure 2: Projections of the two-dimensional unbinned fit to DT events from data onto Mbc (left) and M (right). 


Table 2: The detection efficiencies and the data yields of the ST and DT events. The efficiencies do not include the intermediate 
branching fractions for tt^ 77, 77 77, tt+tt", rj' 7r^7r~r] and rj' —>• p^'y. All uncertainties are statistical only. 


ST mode a 


Kst 


yoT 

K-K- 

47.89 ± 0.35 

1088 ± 40 

13.75 ± 0.14 


K+K-tt- 

44.16 ± 0.18 

5355 ± 118 

12.46 ± 0.14 


K'^K~7r~7r^ 

13.25 ± 0.22 

1972 ± 145 

4.32 ± 0.08 


K^K'^n~7r~ 

24.27 ± 0.37 

595 ± 50 

6.05 ± 0.09 


+ - - 
TT^TT TT 

60.26 ± 0.90 

1657 ± 143 

17.18 ± 0.16 

68±14 

7T~r} 

48.39 ± 0.70 

843 ± 54 

14.82 ± 0.16 


'nTVTVT] 

29.48 ± 0.52 

461 ± 41 

7.91 ± 0.11 



43.11 ± 0.88 

1424 ± 147 

11.96 ± 0.13 


7r~n^r} 

26.02 ± 0.32 

2260 ± 156 

7.90 ± 0.11 



left free in the fit. The ARGUS function of Mbc (a) helps to constrain the description of in 

component ( 6 ). This treatment on background contributions has been verified in MC simulations. There is 
no obvious correlation between Mbc (a) and so the probability density functions (PDFs) of these 

two variables are directly multiplied. We obtain the combined DT yield ?/dt from the unbinned fit shown 
in Fig. 2. Table 2 gives the total yields of DT in data and the corresponding DT efficiencies. Combining 
the yields and efficiencies, we obtain B{Df —^ = ( 8-8 ± 1.8)% with Eq. 4. 


3.2. Measurement of —>• r]'p'^) 

In order to improve the statistical precision, we determine the branching fraction for Df —^ p'p'^ using 
STs. As a standalone measurement, this does not benefit from cancellation of systematic uncertainties as in 
the double-tag method. However, a similar cancellation can be achieved by measuring the signal relative to a 
similar, already well-measured final state. Thus, we measure S(D+ —^ Ty'p'*’) relative to 
using 


B{D+ ^ p'p+)B^^y^BPP^ 

B{Dt ^ K+K-Tr+) ~ 2 / 17 ^"--" ‘ epf 


(5) 


where = B{p'^ —>■ 7 r+ 7 r°)H( 7 r° —>• 77 ). 

The decay Df —>■ K '^is reconstructed in the same manner as reported above in the ST mode. Our 
MC simulation of this mode includes a full treatment of interfering resonances in the Dalitz plot [21]. The 
decay Df —>• rj'p'^ is reconstructed via the decays 77' —>■ 7 r+ 7 r “77 and —>■ tt+tt^, where 7 (tt^) — 77 . We 
apply the same criteria to find 7 r° and 77 candidates as were used in the analysis of Df —>• 77' A. We do not 
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Figure 3: Comparison of the M('7r+7r'^) (left) and (right) distributions in ST events of Dj" —f in data (points) 

and inclusive MC (solid line). The arrows show the signal region. 


require PID criteria on the charged tracks, but instead assume them all to be pions. In the reconstruction 
of p"*" and 77 ', the 7 r+ are randomly assigned. The invariant mass, of the candidate is required 

to be within ±0.170 GeV/c^ of the nominal p+ mass, and the invariant mass of the 77 ' candidate, 
is required to lie in the interval (0.943, 0.973) GeV/c^. Additionally requiring 1.955 < Mbc < 1-985 GeV/cP 
to enrich signal events, the M( 7 r+ 7 r°) distribution of Df —>■ 77 '^+ in inclusive MG simulations and data in 
Fig. 3 show good agreement. The small difference visible in the distribution will be taken into 

account in the systematic uncertainties. 

If multiple r ]'candidates are found in an event, only the one with the smallest | AAl is kept. We require 
—0.035 < AA < 0.023 GeV for data and —0.037 < /S.E < 0.029 GeV for MC. Fits to the Mbc distributions 
are used to extract signal yields. To separate the three body process —>■ ? 7 ' 7 r+ 7 r° from the two body decay 
D'^ —>■ rj'p'^, the helicity angle 6^^+ is used to extract the p'^ component, where 9^^+ is the angle between the 
momentum of the 7 r+ from the p'^ decay and the direction opposite to the D'^ momentum in the rest 
frame. The signal ^ 77 '^+ is distributed as cos^ , while the three body process is flat in cos 9^+. 

We perform a two dimensional unbinned maximum likelihood fit to the distribution of Mbc versus 
cos 07 r+ to determine the yield . The signal model of Mbc is the same as that in the analysis of 
^ 77 'A. For cos0,r+j the signal shapes of —>■ 77 '^+ and ^ ? 7 ' 7 r+ 7 r° are determined based on MC 

simulations. Background contributions in Mbc are modeled with an ARGUS function, while background 
contributions in cos0^+ are taken from the events in the Mbc sidebands 1.932 < Mbc < 1-950 GeV/c^ and 
1.988 < Mbc < 1-997GeV/c^. There is no obvious correlation between Mbc and cos9^+, so the PDFs 
used for these two variables are directly multiplied. Figure 4 shows the projections of the two-dimensional 
fit results in data. In the right plot, we further require 1.955 < Mbc < 1-985 GeV/c^ to enrich signal 
events. The fit returns 7 /gr^ = 210 ± 50, and 77 gr |7 = —13 ± 56, which indicates that no significant 

non-resonant D+ —>■ signal is observed. An upper limit of B{Df —>■ 77 ' 7 r+ 7 r°) at the 90% confidence 

level is evaluated to be 5.1%, after a probability scan based on 2000 separate toy MC simulations, taking 
into account both the statistical and systematic uncertainties. As shown in Fig. 5, we see obvious Df signals 
in the Mbc distribution with the requirement of jcos^jr+l > 0-5, while it is not the case when requiring 
I cos 07 r+l < 0.5. This indicates that the three body process is not significant. 

We study the Mbc distributions for events in p+ and 77 ' sidebands. The p+ sideband region is chosen 
as M( 7 r+ 7 r°) < 0.500 GeV/c^, and the 77 ' sidebands are 0.915 < < 0.925 GeV/c^ and 0.990 < 

< 1.000 GeV/c^. No D+ signal is visible in the sideband events, further substantiating that the 
non-resonant processes D'^ —>■ and D'^ —> p'^ are negligible. A simulation study shows that 

the potential background contribution from 77 ' —^ is negligible. 
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Figure 4: Projection plots of the two dimensional unbinned fit onto Mbc (left) and cos0^-f (right). The signal events are 
enriched by requiring 1.955 < Mbc < 1.985 GeV/c^ in the right plot. 
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Figure 5: Mbc distributions with the requirement of |cos0^+| < 0.5 (left) or | cos0^+1 > 0.5 (right). 
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Table 3: Summary of relative systematic uncertainties in percent. The total uncertainty is taken as the sum in quadrature of 
the individual contributions._ 


Source 

B(Dt ^ JX) 

B{Dt ^V'p+) 

MDC track reconstruction 

2.0 


PID 

2.0 

3.0 

TT^ detection 


2.4 

77 detection 

2.7 

3.5 

HE requirement 

1.0 

1.4 

requirement 
M{ri'TTin]) backgrounds 

1.5 

2.0 

Peaking backgrounds in ST 

0.3 


Mbc signal shape 

1.0 

0.6 

Mbc fit range 

1.7 

0.5 

cos 6^+ backgrounds 


2.9 

Uncertainty of efficiency 

1.6 

0.5 

Quoted branching fractions 

1.7 

3.8 

Total 

5.3 

7.5 


The detection efficiency is estimated to be (9.80 ± 0.04)%. Combined with the results for the 

normalization mode K~^, as given in Table 2, we obtain from Eq. (5) the ratio of B{D'^ v'P^) 
relative to B{Df —>■ K~Tr~^) as 1.04±0.25. Taking the most precise measurement of B{D'^ —/C+ZC“7r+) 
= (5.55±0.19)% from CLEO [11] as input, we obtain B{Df —>• rj'p'^) = (5.8±1.4)%. 

3.3. Systematic uncertainties 

In the measurement of B{D^ —>■ rj'X), many uncertainties on the ST side mostly cancel in the efficiency 
ratios in Eq. (4). Similarly, for Df —>■ r]'^ the uncertainty in the tracking efficiency cancels to a negligible 
level by taking the ratio to the normalization mode 11+ —>■ in Eq. (5). The following items, 

summarized in Table 3, are taken into account as sources of systematic uncertainty. 

a. MDC track reconstruction ejjiciency. The track reconstruction efficiency is studied using a control sample 
of Z1+ —?► itr“7r+7r+ in the data sample taken at y/s = 3.773 GeV. The difference in the track reconstruction 
efficiencies between data and MC is found to be 1.0% per charged pion and kaon. Therefore, 2.0% is 
taken as the systematic uncertainty of the MDC track reconstruction efficiency for £>+ —>■ rj'X. 

b. PID ejjiciency. We study the PID efficiencies using the same control sample as in the track reconstruction 
efficiency study. The difference in PID efficiencies between data and MC is determined to be 1.0% per 
charged pion or kaon. Hence, 2.0% (3.0%) is taken as the systematic uncertainty of the PID efficiency 
for Dj rj'X (D+ ^ ri'p+). 

c. 7r° and rj detection. The 7r° reconstruction efficiency, including the photon detection efficiency, is studied 
using a control sample of D° —ZsT^tt+tt^ in the data sample taken at y/s = 3.773 GeV. After weighting 
the systematic uncertainty in the momentum spectra of tt*^, 2 .8% is taken as the systematic uncertainty 
for the TT^ efficiency in Dj —>■ p'p'^. Similarly, the systematic uncertainty for the p efficiency in Dt^ 
p'X {Dj —>• p'p'^) is determined to be 2.7% (3.5%) by assuming data-MC differences have the same 
momentum-dependent values as for detection. The systematic uncertainties were set conservatively 
using the central value of the data-MC disagreements plus 1.0 (1.64) standard deviations for 7r° (ry), as 
appropriate for a 68% (95%) confidence level. Here we inflate the p uncertainty, because the uncertainty 
of the p detection is estimated referring to 7r°. 

d. AE requirement. Differences in detector resolutions between data and MC may lead to a difference 
in the efficiencies of the AE requirements. In our standard analysis procedure, we apply different AE 
requirements on data and MC, to reduce the systematic uncertainties. To be conservative, we examine 
the relative changes of the efficiencies by using the same AE requirements for MC as for data. We assign 
these changes, 1.0% for Dj —>■ p'X and 1.4% for Dj —>■ p'p'^^ as the systematic uncertainties on the AE 
requirement. 


10 








e. requirement. In the right plot in Fig. 3, the resolution of the rj' peak in MC is narrower than 

data. We take the change in efficiency of 2.0%, after using a Gaussian function to compensate for this 
resolution difference, as the systematic uncertainty of the requirement for Df —>■ rj'p'^. 

f. baekground contributions. In the measurement of B(D'^ —> ^ two-dimensional fit is 

performed to the Mbc(ST) and distributions. The uncertainty due to the description of the 

background contributions is estimated by repeating the fit with higher order polynomial func¬ 
tions. We take the maximum relative change of 1.5% in the signal yields as the systematic uncertainty 
on background contributions. 

g. Peaking background contributions in ST. For the ST D~ candidates, we study the potential peaking 
background contributions with the inclusive MC sample. We find that there is no peaking background 
contributions except for Dj —> 7 r+ 7 r“ 7 r~. We consider the rate of peaking background contributions in 
the ST yields, and take 0.3% as the systematic uncertainty of peaking background contributions in the 
ST events. 

h. Mbc signal shape. To estimate the uncertainty in the Mbc signal shape, we perform alternative 
fits with MC-determined signal shapes with different requirements on the truth matches. We take the 
resultant changes of 1.0% and 0.6% in B{Df v'X) and B{Df —^ v'p^) as the systematic uncertainties, 
respectively. 

i. Mbc fit range. We change the fit ranges of Mbc for ST modes, and take the resulting changes of 1.7% 

and 0.5% in B{Dfi —>• r]'X) and B{Df —>• as the systematic uncertainties, respectively. 

j. cos 07 r+ background contributions. In the measurement of B{Dfi —> rj'p'^), a two dimensional fit is 
performed to the Mbc and cos9t^+ distributions. The shape of the backgrounds in cos 07 r+ is taken from 
the kernel-estimated distribution of the events in the Mbc sidebands with the kernel width parameter 
p = 2 [22] . The uncertainty due to the description of the cos 9.,^+ background contributions is estimated by 
repeating the fit with p = 1.5. We take the relative change of 2.9% in the signal yields as the systematic 
uncertainty on cos 9.„+ background contributions. 

k. Uncertainty of efficiency. In the measurement of B{Dfi —>■ rj'X)^ we use the inclusive MC samples to 

determine e^rp. The DT efficiency erQrp is determined by erQrp = where is obtained 

from MC simulated events of Dj —>■ a and Dfi —>■ 77 '/?, and /3 refers to the five most dominant final states 
7 r+, i4r+, P+, e+z/g and and Bp is the decay rate of Dfi rj'l3. We assign the world averages to the 

branching fractions of these five modes, except for B{Dfi —>■ rj'p'^)^ which is taken from our measurement. 
The statistical uncertainties in sind the Bp uncertainties are propagated to The uncertainties 
of Egx and e^x s-re propagated to B{Dfi —>■ rj'X) and yield a systematic uncertainty of 1.6%. For the 
measurement of B{Dfi —>■ 'ij'p'^), the uncertainty of the efficiency due to the limited MC statistics is 
estimated to be 0.5%. 

l. Quoted branching fractions. The branching fractions of rj' —>■ mrrj, rj —>■ 77 , —>■ 77 are taken from 

PDG [9]; the branching fraction for Dfi —>■ K'^K~7t~^ is taken from CLEO’s measurement [11]. Their 
uncertainties are 1.6%, 0.5%, 0.03% and 3.4%, respectively. 


4. Summary and Discussion 

We measure the branching fraction B{Dfi —>■ rj'X) = ( 8.8 ± 1.8 ±0.5)%, which is consistent with CLEO’s 
measurement [10]. The weighted average of these two results is B{Df —>• rj'X) = (10.3 ± 1.3)%. We also 
measure the ratio B{Dfi —>■ rj'p'^) / B{Dfi —>■ iL+iL“ 7 r+) = 1.04 ± 0.25 ± 0.07, from which we get B{Df —>■ 
? 7 'p+) = (5.8 ± 1.4 ± 0.4)%. This is nearly half of CLEO’s older result [7], but compatible with CLEO’s 
newer measurement of B{Dfi —>■ p'tt+tt*^) [11], in which the resonant process rj'p'^ is believed to dominate. 
We also report a limit on the non-resonant branching ratio B{Dfi —>■ p'tt+tt^) < 5.1% at the 90% confidence 
level. These results reconcile the tension between experimental data and theoretical calculation [8]. Taking 
the world average values of other exclusive branching fractions involving rj' as input, we obtain the sum of 
exclusive branching fractions B[Dfi —>• rj'K'^, ij'n'^,rj'p'^,rj'lvi) = (11.9± 1.6)%, in which I denotes e+ or p+, 
and where we have assumed that B{Df —>• rj'= B{Df —^ p'e+z/g)- This summed exclusive branching 
fraction is compatible with the new weighted inclusive result B{Df —>• ij'X) = (10.3 ± 1.3)%. 
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